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Abstract 

We derive model equations for optical pulse propagation in a medium 
described by a two-level Hamiltonian, without the use of the slowly varying 
envelope approximation. Assuming that the resonance frequency of the two- 
level atoms is either well above or well below the inverse of the characteristic 
duration of the pulse, we reduce the propagation problem to a modified 
Korteweg-de Vries or a sine-Gordon equation. We exhibit analytical solutions 
of these equations which are rather close in shape and spectrum to pulses in 
the two-cycle regime produced experimentally, which shows that soliton-type 
propagation of the latter can be envisaged. 



1 Introduction 



Recent advances in dispersion managing now allow the generation of ultra- 
short optical pulses containing few oscillations directly from a laser source. 
Two-cycle pulses have been recently reported in mode-locked Ti-sapphire 
lasers using double- chirped mirrors [U El IS] ■ Because the pulse duration be- 
comes close to the optical period, a question of interest is to know if few cycle 
pulses can generate optical solitons in nonlinear media. The usual description 
of short pulses propagation in nonlinear optics is made using the nonlinear 
Schrodinger (NLS) equation which is derived using the slowly varying enve- 
lope approximation. However, for ultrashort pulses considered in this paper 
the slowly varying envelope approximation is not valid any more. This situa- 
tion, and the corresponding one in the spatial domain, called 'non-paraxial', 
have already given rise to several studies. An approach consists in adding 
corrective terms to the NLS model [3]. This high-order perturbation ap- 
proach still involves the slowly varying envelope approximation, and requires 
cumbersome and difficult computations. The approach of [5] allows one to 
determine the ray trajectories in a very rigorous way, without any use of the 
paraxial approximation. However, it still makes use of the slowly varying 
envelope approximation in the time domain, and therefore can hardly be 
generalized to the problem under consideration in this paper. It is prefer- 
able to leave completely the concept of envelope. Indeed, it is not adapted 
when a pulse is composed of few optical cycles. The aim of this paper is to 
demonstrate that other approximations can be envisaged and can also lead 
to completely integrable equations, and support solitons. The basic principle 
of our work is that a soliton can propagate only when the absorption is weak, 
therefore its characteristic frequency must be far away from the absorption 
range of the material. If it is far below, a long-wave approximation can be 
performed. On the other hand, if it is far above, it will be a short-wave ap- 
proximation. Both approaches are used in this paper leading to completely 
integrable systems. It is organized as follows. In section 2 we develop the 
model which is based on a non absorbing homogeneous and isotropic two- 
level medium. A semi-classical approach is used leading to the well-known 
Maxwell-Bloch equations. The long-wave approximation is investigated in 
section 3. In this case the model reduces to a modified Korteweg-de Vries 
(mKdV) equation. The two-soliton solution is very close to the experimen- 
tally observed two-cycle pulses. In section 4 we investigate the short-wave 
approximation. The resulting model is formally equivalent to that describing 
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the self- induced transparency. It can be reduced to the sine- Gordon equa- 
tion. Again, the two-soliton solution is comparable with the experimental 
observations 

2 Model 

In this section we derive the starting equations for further analysis. The 
medium is treated using the density matrix formalism and the field using the 
Maxwell equations. 

We consider an homogeneous medium, in which the dynamics of each 
atom is described by a two-level Hamiltonian 

A more realistic description should take into account an arbitrary number of 
atomic levels. Indeed, we consider wave frequencies far from the resonance 
line of the medium, and in this situation all transitions should be taken into 
account. But we intend here to suggest a new approach to the description 
of ultrashort optical pulses. Therefore we restrict the study to a very simple 
and rather academic model. 

The atomic dipolar electric momentum is assumed to be along the x-axis. 
It is thus described by the operator jl = pe x , where e x is the unitary vector 
along the x-axis and 

"=(;s)- < 2 > 

The polarization density P is related to the density matrix p through 

P = Ntr(pii), (3) 

where N is the number of atoms per unit volume. Thus P reduces to Pe x . 

The electric field E is governed by the Maxwell equations. In the absence 
of magnetic effects, and assuming that the wave is a plane wave propagating 
along the z-axis, polarized along the x-axis, E = Ee x , they reduce to 

d 2 z E=\d?(E + 4nP). (4) 

c is the light velocity in vacuum. We denote by d t the derivative operator 
with regard to the time variable t, and so on. 
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The coupling between the atoms and the electric field is taken into account 
by a coupling energy term in the total Hamiltonian H, that reads: 



H = H - fiE. (5) 
The density matrix evolution equation (Schrodinger equation) writes as 

thd tP =[H,p] + n, (6) 

where 1Z is some phenomenological relaxation term. The set of equations 
(JIUO]) is sometimes called the Maxwell-Bloch equations, although this name 
denotes more often a reduction of it. 
Setting 

t' = c t , P' = AnP , p 1 = ArrNhcp , a' = (7) 

he 

it I H Til _ H J _ U a ,b ( o\ 

allows one to replace the constants c, N, h and 47r in system (JBHOJ) by 1. We 
denote the components of p by 

'-(as)- « 9 » 

and so on, and by O = u>b — uj a the resonance frequency of the atom. 
The relaxation expresses as 

n = ih( phl Jt ~ pt/ , Tt ), (io) 



-Pt/n -Pb/n 

where Tj, and r t are the relaxation times for the populations and for the coher- 
ences respectively. We show below that, according to the fact that relaxation 
occurs very slowly with regard to optical oscillations, the relaxation term 1Z 
could be omitted. 



3 Long-wave approximation 

3.1 A modified Korteweg-de Vries equation 

Let us first consider the situation where the wave duration t w is long with 
regard to the period t r = 2n/Q (recall that Q = uj^ — u a ) that corresponds to 
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the resonance frequency of the two-level atoms. We assume that t w is about 
one optical period, say about one femtosecond. Thus we assume that the 
resonance frequency Q is large with regard to optical frequencies. In order to 
obtain soliton-type propagation, nonlinearity must balance dispersion, thus 
the two effects must arise simultaneously in the propagation. This involves 
a small amplitude approximation. Further, we can speak of soliton only if 
the pulse shape is kept on a large propagation distance. Therefore we use 
the reductive perturbation method as defined in j^j. We expand the electric 
field E, the polarization density P and the density matrix p as power series 
of a small parameter e as 

e=J2 £HE - > p = Y. £np - * p = 5>>«' ( n ) 

njjl n^l riJSO 

and introduce the slow variables 

T = e(t-y) , C = e 3 z. (12) 

Expansion (|TT| gives an account of the small amplitude approximation. The 
retarded time variable r describes the pulse shape, propagating at speed 
V in a first approximation. Its order of magnitude e gives account for the 
long-wave approximation, so that the pulse duration t w has the same order 
of magnitude as t r /e. The propagation distance is assumed to be very long 
with regard to the pulse length ct w , therefore it will have the same order 
of magnitude as ct r /e n , where n ^ 2. The value of n is determined by the 
distance at which dispersion effects occur. According to the general theory 
of the derivation of KdV-type equations jSj, it is n = 3. The ( variable of 
order e 3 describes thus long-distance propagation. The physical values of the 
relaxation times T\, and r t are in the picosecond range, or even slower, thus 
very large with regard to the pulse duration t w . Therefore we write 

Tj = -4: for j = b and t. (13) 

e l 

The Schrodinger equation © at order e° is satisfied by the following 
value of po, which represents a steady state in which all atoms are in their 
fundamental state a: 

0»=(o°o)- (") 
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Notice that, according to the change of variables (J7J), the trace tr(p) of the 
density matrix is not 1 but a = AirNhc. Then the Schrodinger equation (jBJ) 
at order e 1 yields 

Pit = (15) 

so that 

Pl = Ma*. (i6, 

The Maxwell equation (jlj) at order e 3 gives the value of the velocity 

in accordance with the limit of the dispersion relation as the frequency uo 
tends to zero. 

The Schrodinger equation © at order e 2 yields p\ a = pxb = and 

= (18) 

Then 

P 2 = (19) 

and the Maxwell equation (jlj) at order e 4 is automatically satisfied. 
The Schrodinger equation (jUJ) at order e 3 gives 

P2b = ~P2a = lj ^E 2 (20) 

and 

P3t " TT 3 _ Tp" T " Yf^ 1 _ _ n^ 1 ' (21) 

The corresponding term of the polarization density P contains a nonlinear 
term: 

Ps - -q~ E3 ~ -qt^ Ei - ~or Ei ' (22) 

but the terms involving the relaxation do not appear. The Maxwell equation 
at order e 5 yields the following evolution equation for the main electric field 
amplitude E\\ 

d(El = + 2 -^S T El (23) 



which is a mKdV equation. Equation (J23|) can be generalized as follows: a 
general derivation of KdV-type models [7j shows that the coefficient of the 
dispersive term d 3 E\ in this equation must be (l/6)d 3 k/du 3 . We check by 
direct computation of the dispersion relation that it holds in the present case. 
Another heuristic reasoning can relate the value of the nonlinear coefficient 
of equation (J23j) to the third order nonlinear susceptibility x ■ It uses the 
nonlinear Schrodinger (NLS) equation which describes the evolution of a 
short pulse envelope in the same medium. The NLS equation writes as ( 8 , 
6.5.32) 

id ( £ - ~k 2 c%£ + 7 £ \Ef = 0, (24) 

where £ is the envelope amplitude of the wave electric field, ki the group 
velocity dispersion, and 7 is related to x^ through 

QUI! ,nN , , 

7 = X (3) - 25 

nc 

n is the refractive index of the medium, u the wave pulsation. We drop the 
dispersion term, replace uo by id T , and notice that £ coincides with the real 
field E\ at the long- wave limit, to get 

dfr = —X^drEl (26) 
nc 

Equation (|26|) gives an expression of the nonlinear coefficient in the mKdV 
equation (|23|) . The relevant component of the third order nonlinear suscep- 
tibility tensor x^ computed from the above model is jH] 

XtttAVi ^1 u i ~ UJ ) — -^77 9 • (27) 

Taking the long- wave limit u — > in equation ()27jl. we check that the 
expression of the nonlinear coefficient obtained from (J26|) holds in the present 
case. Thus we can write equation (|2*3*|l as 



d ' El = 6 J? 



drEi Xxlxxiu, w, u, -u) 



d T E\ . (28) 

uj=0 



It can be reasonably conjectured that equation (j28j) will still hold in the 
more general case of an arbitrary number of atomic levels, when the inverse 
of the characteristic pulse duration is much smaller than any of the transition 
frequencies of the atoms. 
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3.2 The two-soliton solution 



The mKdV equation ()23j) is completely integrable by means of the inverse 
scattering transform jU]. The iV-soliton solution has been given by Hirota 
[TUj . In order to write it easily, we write the mKdV equation (}2"3*)) into the 
dimensionless form 

d z u + 2d T u 3 + <9|u = 0, (29) 

where u is a dimensionless electric field, and Z and T dimensionless space 
and time variables defined by 

u=lr , Z=^- T = ^-. (30) 
The characteristic electric field, space and time are defined by 

in normalized units. Relations (|3UH31|) can be expressed in a more convenient 
way as follows. Let us first choose as reference time the pulse length t w (in 
physical unit). The small perturbative parameter is then 

e = sir <32) 

The characteristic electric field S, and propagation distance C, are 

h hc 2 QH 3 

S = tM ' £ = w#' (33) 

where the speed V is 

Then the quantities involved by the dimensionless equation (J29|) are related 
to the quantities measured in the laboratory through 

-=! • z =t • r 4(<-^)- < 35 > 
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The soliton solution writes as 



u = psechrj, 



(36) 



with 



•q = pT — p 3 Z — 7, 



p and 7 being arbitrary parameters. 
The two-soliton solution is 



(37) 



+ e" 2 + 



yi-yi \ 

P1+P2 J 



1 + 



+ 



4p? (P1+P2) 



e 2 "2 
4pJ 



+ 



Pi— P2 
P1+P2 



4 ' 
e 277 1 +2i 72 

16pfpl 



with 



m = pj t - Pj z - 7i 



(38) 



(39) 



for j = 1 and 2. The parameters pi, p2> 7i an d 72 are arbitrary. When 
they take real values, the explicit solution (|58*|l describes the interaction of 
two localized bell-shaped pulses, which are solitons. An example of this 
solution is drawn on figure [TJ using the values of the parameters p% = 3, 
P2 = 5, 71 = 72 = 0. But expression (jHEj) also describes the so-called higher 




Figure 1: Two-soliton solution of the mKdV equation, using dimensionless parameters. 
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order solitons, which can be considered as a pair of solitons of the above 
kind linked together, and have often an oscillatory behaviour. An example 
is given on figure 121 It uses the values of parameters p\ = 1 + 4i, p 2 = pi, 
71 = —72 = in/2. The corresponding spectrum is drawn on figure El 




Figure 2: Second-order soliton solution of the mKdV equation, using dimensionless 
parameters. 

These spectrum and pulse profile are comparable to the experimental pulses 
given by [3]. It can thus be thought that the two-cycle pulses produced 
experimentally could propagate as solitons in certain media, according to 
the mKdV model. 

4 Short-wave approximation 
4.1 A sine-Gordon equation 

We now consider the situation in which the resonance frequency Q of the 
atoms is below the optical frequencies. Then the characteristic pulse dura- 
tion t w is very small with regard to t r = 2it/VL, thus we use a short-wave 
approximation. We introduce a small perturbative parameter e, such that 
the resonance period t r = t r /e, where t r has the same order of magnitude 
as the pulse duration t w . The perturbative parameter e is thus about t w /t r . 
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Consequently, the Hamiltonian H$ of the atom is replaced in the Schrodinger 
equation (0) by 

eH . (40) 
We introduce a retarded time r and a slow propagation variable ( such that 

r=(f-£) , C = «- (41) 

The zero order reference time is chosen to be t w , therefore r is not a slow 
variable. The definition of the variable ( gives account for long distance 
propagation. Computation shows that the dispersion effects arise at distances 
about ct w /e, from which follows the choice of the order of magnitude of (. 
The electric field E is expanded as E = ^ n> $£ n E n) and so on. The pulse 
duration t w is still assumed to be about one femtosecond, corresponding to 
an optical pulse of a few cycles. The relaxation times 77,, r t are very long 
with regard to t w . Since the above scaling uses t w as zero-order reference 
time, this can be expressed by setting 

7-. 

Tj — — for j = b and t. (42) 
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Notice that (|42j) differs formally from the assumption (J13|) written in the 
previous section but represents the same physical hypothesis. 

The above scaling can also be presented from another viewpoint, taking 
the characteristic time of the resonance as zero-order reference time, as 
follows. The relevant component of the third order nonlinear susceptibility 
tensor x computed from the above model is given by formula (J27j) . where lo 
is the wave frequency (while Q = ujf,—u a is the resonance frequency of the two- 
level system). The short-wave approximation corresponds to uo — > oo. Then 
Xxxxx tends to zero. Thus a linear behaviour of the wave can be expected in 
the short-wave approximation, except if the nonlinearity is very strong. The 
latter physical assumption can formally be expressed by assuming that the 
product pE is very large with regard to Q, according to 

H = H + ~pE, (43) 

where the small parameter e tends to zero. Then the short wave approxima- 
tion can be sought using the expansions 

3^0 j^O j^O 

and slow variables ( and r such that 

d t = -d T , d z = ^-d T + d ( . (45) 
e eV 

The definition of variables ()45|) is very close to the standard short wave 
approximation formalism developed e.g. in JU^]. It is easily checked that 
the scalings (j4UH41j) and (|43H45j) are equivalent. We refer to the former below. 
The Schrodinger equation (jUJ) at order e° yields 

id T p 0a = -E (/j,p* ot - p* p ot ) , (46) 
idrpob = +E (pp* ot - p*p ot ) , (47) 
id T pot = -Eo/i (p ob - p 0a ) . (48) 

From (J46H47)) we retrieve the normalization condition of the density matrix 
<9 T tr p = 0. We introduce the population inversion w = pob — Poa and get 

. (a-w)/2 l pfE w 
P ° 1 -ifi*f T E w (a + w)/2 1 
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(as above, tip = a = 4nNhc due to the normalization), and the equation 

d T w = -%| 2 £ J E w. (50) 

Then expression (J3J) of the polarization P yields Pq = 0, and the Maxwell 
equation (J3J) at order e° becomes trivial if the velocity is chosen as V = 1. 
The Schrodinger equation (JHJ) at order e writes then as 

id rPl = [H , po] - frEo, Pl ] - \fiEk, po)+i( P ° b J T , b ; ~ P ^Jl ) . (51) 

V ~Pot/ T t -pob/n j 

Defining W\ = pu — pi a , the off-diagonal components of equation (|5Tj) yield 
Pi t = i(n + ^jJ Pat + W J (E w 1 + E l w) ) (52) 
so that the corresponding term Pi = p*pu + pp\ t of the polarization is 

Pi = -2VL \p\ 2 / / Ew. (53) 



Notice again that the relaxation does not appear in the expression of the 
polarization at this order. The Maxwell equation (@J) at order e then reduces 
to 

d ( d T E = n\p\ 2 E w . (54) 
Equations ()50I54|) yield the sought system. If we set 



p=-i\p\ / E w, (55) 



they reduce to 



d ( E = ittp, (56) 
d T p = — i\p\ 2 E w, (57) 
d T w = -AiE p, (58) 

which coincide with the equations of the self-induced transparency, although 
the physical situation is quite different: the characteristic frequency l/t w 
of the pulse is far above the resonance frequency Q, while the self-induced 
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transparency occurs when the optical field oscillates at the frequency f2. The 
quantities E and w describe here the electric field and population inversion 
themselves, and not amplitudes modulating a carrier with frequency Q. No- 
tice that E and w are here real quantities, and not complex ones as in the 
case of the self-induced transparency. Further, p is not the polatization den- 
sity, but is proportional to its r-derivative. Another difference is the absence 
of a factor 1/2 in the right-hand side of equation (J56j) . 

Since they explicitly involve the population inversion, the model equa- 
tions (J56H58j) cannot be generalized easily to more realistic situation in which 
an arbitrary number of atomic levels are taken into account. Recall that, ac- 
cording to the assumption made at the beginning of the section, this model 
is valid for a very strong nonlinearity only. In particular, we assumed that 
the atomic dipolar momentum fi has a very large value. In a more realistic 
situation, it can be expected that only the transition corresponding to the 
largest value of the dipolar momentum will have a significant contribution. 
If several transitions correspond to large values of the dipolar momentum 
with the same order of magnitude, we can expect that the short-wave ap- 
proximation will yield some more complicated asymptotic system involving 
the populations of each level concerned. The derivation of such a model is 
left for further study. 



4.2 The two-soliton solution 

Using dimensionless variables defined by 

e = f , w = ± , r = f , z = f, (59) 

E r w r T L 

where the reference values satisfy 

E r To\p\ = 1 and w r LQT \n\ 2 = 2, (60) 

and setting r\ — J Z W, the system (j50H54j) reduces to 

d z d T Q = 2Qd z r], (61) 
d z d T V = -2&d z @. (62) 

Equations (j61H62J) have been found to describe short electromagnetic wave 
propagation in ferrites, using the same kind of short-wave approximation ^T] . 
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The following change of dependant variables: 



dzf] = Acosu, (63) 
d z B = Asinu, (64) 

transforms equations (j61H62|) into [TTJ [T3] 

d T A = 0, (65) 
dzdxu = 2Asinu. (66) 

Since, according to ([63)1. A is a constant, (pj)) is the sine-Gordon equation. 
Before we recall some properties of the latter, let us determine the physical 
meaning of the constant A in the present physical frame. Using relations 
(|63H64|) and the definition of 77, we find that 

A 2 = lim (W 2 + (d z e) 2 ) (67) 

Since is the dimensionless wave electric field, it vanishes at infinity. Thus 
we can have a non- vanishing solution only if some initial population inversion 
Wi = WiW r is present. The constant involved by equation ()66|) is then A = 
W t . 

Using the variable Z = 2WiZ, equation ((HBJ reduces to the sine-Gordon 
equation 

d^dxu = sinw. (68) 

The quantities involved by equation (|55|l are related to the quantities mea- 
sured in the laboratory through 

M • T =H*-D- (69) 



z 



E=—J sinw , w = WiCosu. (70) 
The electric field and propagation length scaling parameters are 

Er = 7^-, (71) 

(72) 



^47rA^|/i| 2 w7i 
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in which the initial population inversion Wi and typical pulse duration t w 
are given. The small perturbative parameter e can be identified with Qt w , 
expressing the fact that t w is very small with regard to 

The sine-Gordon equation ()68|) is completely integrable A iV-soliton 
solution can be found using either the 1ST or the Hirota method. As in 
section 2, we will consider here the two-soliton solution only, which is [T3] 

u = 2iln(4-J, (73) 



with 



./ 



/ = 1 + ie^ + ie^ - |^|C e 7?1+??2 , (74) 



where 

Vj = k j T + -+ lj for j = 1,2, (75) 

ki, 71 and 72 being arbitrary parameters. When they take real values, for- 
mulas (|7HH75j) describe the interaction of two solitons. The behaviour is very 
close to that of the typical two-soliton solution of the mKdV equation shown 
in figure ^ As in the case of the long wave approximation, the two-soliton 
solution (|73H75|) is also able to describe soliton-type propagation of a pulse in 
the two-cycle regime. The corresponding analytic solution is a second-order 
soliton or breather, which can be considered as two bounded solitons, and is 
obtained using complex conjugate values of the soliton parameters k\ and k 2 . 
An example is given in figure EJ with the values of parameters k± = 1 + 4i, 
hi — 1 — 4i, 71 = 72 = 0. The pulse profile, with the corresponding popula- 
tion inversion and spectrum are drawn in figure El The profile and spectrum 
are comparable with the experimental observation of |3 3 . Notice again that 
an initial population inversion Wi 7^ is required. Total inversion (wi = 1) is 
not necessary but, as shows the expression (|72*|) of the propagation reference 
length L, a small inversion reduces the soliton amplitude and increases the 
propagation distance at which nonlinear effects occur. 



5 Conclusion 

We have given two models that allow the description of ultrashort optical 
pulses propagation in a medium described by a two-level Hamiltonian, when 
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Figure 4: Electric field evolution of the second-order soliton solution of the sine-Gordon 
equation, using dimensionless parameters. 



the slowly varying envelope approximation cannot be used. Using approxima- 
tions based on the hypothesis that the resonance frequency of the medium is 
far from the field frequency, we derived completely integrable models. When 
the resonance frequency is well above the inverse of the typical pulse width of 
about one femtosecond, a long-wave approximation leads to a mKdV equa- 
tion. When in the contrary the resonance frequency is well below the field 
frequency, a short-wave approximation leads to a model formally identical 
to that describing self-induced transparency, but in very different validity 
conditions. It can be reduced to the sine-Gordon equation. The scaling 
parameters for these approximations have been written down explicitly. 

Both the mKdV and the sine-Gordon equations are completely integrable 
by means of the 1ST method and admit iV-soliton solutions. The two-soliton 
solution is able to describe the propagation of a pulse in the two-cycle regime, 
very close in shape and spectrum to the pulses of this type produced exper- 
imentally. It does not mean that the formulas of this paper describe the 
experimental results, because we have considered a propagation problem, 
and experimental results concern pulses generated directly at the laser out- 
put. But we have shown that soliton-type propagation, with only periodic 
deformation of the pulse during the propagation, may occur for such type of 
pulses, under adequate conditions. In the short-wave approximation, these 
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Figure 5: (a) Pulse profile, (b) population inversion and (c) spectrum, of the second-order 
soliton solution of the sine-Gordon equation of figure Dimensionless parameters. 



conditions involve an initial population inversion, at least a partial one. 

Further, the study of a two-level Hamiltonian can be considered as an 
academic problem, showing the tractability of such an approach. A rather 
remarkable feature is that the computations involved by the derivation of 
the asymptotic models are relatively short and easy. Therefore the appli- 
cation of the same approach to more realistic situations can be reasonably 
envisaged. A generalization of the mKdV equation obtained in the long-wave 
approximation has been proposed on heuristic grounds, and should be jus- 
tified rigorously. A generalization of the model obtained in the short-wave 
approximation would require a special study. Last, in a more realistic model, 
it can be envisaged that some transition frequencies are well above the inverse 



17 



of the characteristic pulse duration, but that some other are below it. The 
treatment of such a situation will mix the above short-wave and long-wave 
approximations, it is left for further study. It can be expected that the result 
will depend strongly on the particular physical situation considered. 
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